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Abstract 

One of the key merits of VT— symmetric (i.e., parity times time reversal 
symmetric) quantum Hamiltonians H lies in the existence of a horizon of the 
stability of the system. Mathematically speaking this horizon is formed by 
the boundary of the domain C 1R D of the (real) coupling strengths for 
which the spectrum of energies is real and non-degenerate, i.e., in principle, 
observable. It is shown here that even in the elementary circular four-site 
quantum lattices with D = 2 or D = 3 the domain of the hidden Hermiticity 
proves multiply connected, i.e., topologically nontrivial. 



1 Introduction 



One of the most interesting formulations of the standard and robust dictum 
of quantum mechanics emerged in connection with the acceptance of the so 
called VT— symmetric operators of observables where V means parity while 
T represents time reversal (cf. review papers [TJ |2j [3] for an exhaustive 
discussion). One of the main reasons of the last-year rebirth of interest in 
this new paradigm may be seen, paradoxically, in its impact on classical 
experimental optics [3]. 

The latter experimental activities (i.e., basically, the emergence of a 
few successful classical-physics simulations of quantum effects) re-attracted 
attention to the innovative theory. We may mention, pars pro toto, pa- 
per j5] which offered an exhaustive constructive classification of all of the 
VT— symmetric quantum Hamiltonians H defined in the finite- dimensional 
Hilbert spaces H of dimensions N = 2 and N = 3. 

Inside the most elementary N < 3 family of models no real surprises 
and spectral irregularities have been encountered. In contrast, in Ref. [6] 
we found that certain anomalies certainly emerge at N = 8. In our present 
brief continuation of these developments we intend to show that the sim- 
plest models exhibiting similar irregularities in their spectra already occur, 
unexpectedly, as early as at the next Hilbert-space dimension N = A. 

2 The four-site quantum-lattice models 
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Figure 1: Graphical symbol for the straight-line open-end four-site lattice. 
The numbers in the small circles (= sites) are the unperturbed energies while 
the letters b and c near the nearest-neighbor-interaction lines represent the 
(real) couplings. 



2.1 The exactly solvable straight-line case. 

In Refs. [7] the successful tractability of more-than-three-dimensional Hamil- 
tonian matrices resulted from a drastic simplification of their structure. We 
merely admitted their tridiagonal versions. In the language of physics this 
corresponded to the picture in which the system lived on a N— site straight- 
line lattice endowed with the mere nearest-neighbor interactions. At iV = 4 
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this is schematically depicted in Fig. [TJ The small circles represent the sites 
while their frame-line connections symbolize the interactions. 

The left-right symmetric straight-line lattice of Fig. [1] (i.e., of Refs. [7]) is 
being assigned the Hamiltonian given in the form of the two-parametric real 
matrix 
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The quantitative analysis of the models of this form is more or less trivial 
even at larger N > 4. Curious reader may find many details, say, in review 
paper [8]. 

2.2 VT— symmetric circular lattices and their simplest 
four-site example 




Figure 2: The circular four-site lattice. 



Once we replace Fig. [T] by its circular version of Fig. [2] we may immediately 
interpret the new diagram as representing the new N = 4 quantum model 
which is given by the following three-parametric four-by-four matrix form of 
the Hamiltonian as proposed and studied in Ref. [9], 
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(2) 



It may be shown that in the new model with one more coupling which con- 
nects the "upper two" sites, the method of construction of the boundary 
does not change. At any number N of sites along the (circular) lattice 
the reality property of the spectrum of the energies will remain tractable by 
the standard mathematical techniques, a few representative samples of which 
may be found in Ref. [10]. Interested readers may also search for a broader 
mathematical context, say, in Refs. [TT| IT2]. 
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Once we restrict our attention just to our special toy model H^ 4 \a, b, c) it 
proves sufficient to recall the entirely elementary considerations of Ref. [TO] . 
This leads to the conclusion that the spectrum of energies stays real and 
nondegenerate if and only if the triplet of parameters (a, b, c) lies inside the 
domain 

D (H) := {(a,b,c) E R 3 \ W(a,b,c) > 0, Q(a,b,c) > 0, P(a,b,c) > } 

(3) 

where 

W(a, b, c) = (8 + c 2 - a 2 ) 2 - 4 [16 - (a + c) 2 ] b 2 , (4) 
Q(a, 6, c) = [(a + 3)(c - 1) - 6 2 ] [(a - 3)(c+ 1) - 6 2 ] (5) 

and 

P(a,b,c) = 10- a 2 - 2 6 2 - c 2 . (6) 

In other words, for the couplings moving to the two-dimensional surfaces of 
J)( H ) from inside one observes that the quadruplets of the real bound-state 
energies themselves behave in an easily understandable manner. The reason 
is that we may rewrite the secular equation in the form S(s, a, b, c) = where 
the energies E± = ±\/s emerge in pairs and where 

S(s, a, b,c) := s 2 + (-10 + c 2 + 2 b 2 + a 2 ) s+9+6 6 2 -9 c 2 +6 4 -2 cab 2 -a 2 +c 2 a 2 . 
This recipe generates the two auxiliary roots 

4s = 4s (±) = 20 - 2 a 2 - 2 c 2 - 4 b 2 ± 2 y/W(a, b, c) 
where we already know the function of Eq. (J3J), 

W(a, 6, c) = 64 + 16 c 2 - 64 6 2 - 16 a 2 + c 4 + 4 c 2 6 2 - 2 c 2 a 2 + 4 b 2 a 2 + a 4 + 8 cafe 2 . 
In the spirit of the general results of Ref. [9] we may summarize that 

1. whenever W(a, b, c) — > + one can spot the two pairs of energies which 
approach the two distinct values 

E (w=o) = ±v /(i -a 2 -2 6 2 -c 2 )/2 
representing the two limiting doubly-degenerate energies; 

2. whenever Q(a,b,c) — > + just the two energies move to zero while the 
other two energies do not vanish in general, 

E ,3 = ±VlO-a 2 -2b 2 -c 2 ; 

3. for P(a, b, c) — > + we must expect that all of the four real energies will 
vanish simultaneously. 
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Figure 3: The graphical determination of the innermost star-shaped domain 
T>( H > as assigned to the quantum lattice of Fig. [TJin Refs. [TJ. 



3 Two-parametric simplified versions of the 
circular four-site lattice 

3.1 The case of a = 

Naturally, in the no-upper-interaction limit lim a _j. (a, b, c) = (b, c) 
we return to the elementary straight-line model of Fig. [TJ For our present 
purposes it is then sufficient to recall that the main features of such a sim- 
plified model were described in Refs. [TJ. In particular, we know that at 
N = 4 the spectrum of energies remains real and nondegenerate inside the 
innermost star-shaped domain T>( H ' c R 2 shown, in Fig. [3j as lying inside 
an auxiliary circumscribed ellipse. The boundary dT>^ H > (i.e., the physical 
horizon of the system in question) is composed of the four hyperbola-shaped 
curves. The key features of this example (like the triple intersections of the 
boundaries, etc.) generalize, mutatis mutandis, to the family of the similar 
models at all of the dimensions N < oo [5]. 

We are now prepared to replace the elementary and transparent graphical 
determination of the star-shaped domain P^' assigned to the straight-line 
quantum lattice and displayed in Fig. [3] by its much more complicated o ^ 
analogue. At a freely variable a the knowledge of the a = section may 
and will still serve us as a very useful independent test of our forthcoming 
observations and conclusions. 

3.2 The case of b = 

In another preparatory step we notice that our present four-site toy model 
H^(a, b, c) degenerates to the trivial non-interacting composition (i.e., direct 
sum) of the two N = 2 models at b = 0. For this reason the 6 = limiting 
case should be considered exceptional. 

It is worth emphasizing that at the 6 = two-dimensional special case, 
even the general definition of the domain T>^ is slightly misleading. Indeed, 
Fig. H] which displays the three sets of boundaries (viz., the two hyperbolas 
W(a, 0, c) = 0, the four straight lines Q(a, 0, c) =0 and the single circle 
P(a, 0, c) = 0, respectively) should not be taken too literally because one of 
the boundaries (viz., the doublet of hyperbolas W(a, 0, c) = 0) describes in 




Figure 4: The degenerate case of the simply connected rectangular domain 
at b = 0. 

fact a sign-non-changing (i.e., the reality-of-energies non-changing) curve of 
the doubly degenerate (and, hence, irrelevant and removable) zeros of the 
function W(a, 0, c) = (8 + c 2 - a 2 ) 2 . 

This means that at b = the domain T)^ is strictly rectangular and 
strictly simply connected. In this context one of the key messages of our 
present study is the surprising discovery of the loss of both of these properties 
in the general case with the freely variable coupling strength b. 

Incidentally, the multinomial W(a, b spec , c) becomes factorizable also at 

bspec 1) 

W(a, 1, c) = (a + c) (a 3 -ca 2 -12a- c 2 a + 20 c + c 3 ) . 

This is an artifact which does not carry any immediate physical meaning. Its 
manifestation is of a purely geometrical character which will only be briefly 
mentioned later. 




Figure 5: The triply connected nature of the triple-overlap domain T>^ H > 
for the quantum lattice of Fig. [2] at c = (i.e., in the no-central-coupling 
extreme). 

3.3 The case of c = 

In the third (and last) preparatory step let us discuss the vanishing-coupling 
special case in which c = and 

W(a, b, 0) = (8 - a 2 ) 2 - 4 [16 - a 2 ] b 2 , (7) 
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Figure 6: Same as Fig. [5] (detail). 




Figure 7: The deformation of Fig. [6] at c = 0.95. 



and 



Q(a, b, 0) = (3 + b 2 + a) (3 + b 2 - a) 



P(a,b,0) = 10 -a 2 -2b 2 



(9) 



The detailed study of precisely this special case reveals in fact the possibility 
of the emergence of a topological nontriviality in the general case. The 
detailed form of such a c = hint may be seen in Fig. |5] and in its magnified 
version [6j Indeed, as long as the condition Q(a,b, 0) > degenerates to the 
elementary constraint 

3 + b 2 > a > -3 - b 2 

just the left and right small horizontal-parabolic segments (with their extreme 
at b = and \a max \ = 3) should be cut out of the elliptic domain with 
P(a, b, 0) > as inadmissible since Q(a, b, 0) < there. As long as we only 
have 16 > a 2 , the remaining constraint W(a, b, 0) > acquires the elementary 
form 



l&i<7 



— a 



4 V16 - a 2 ^ 

of the geometric limitation of the admissible range of b by the two intersecting 
or rather broken and touching curves. Indeed, one must keep in mind that 



the nonnegative function W(a, 0, 0) = (8 



,2\2 



solely vanishes at a — 8. 



These observations imply that the allowed region decays into the three dis- 
connected open sets (cf. Figs. |5]-[7j). The big one is formed by the eye-shaped 
vicinity of the origin, with its extremes at the points (a, b)± = (±v8, 0). The 
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other two smaller open sets are fish-tail-shaped. In the pictures (cf., in partic- 
ular, the magnified right one in Fig. [6]) these two domains are easily spotted 
as containing the respective 6 = intervals of \a\ E (y/8, 3). In this sense they 
may be expected to support a perturbatively inaccessible "strong-coupling" 
dynamical regime. 

An additional indication of the suspected emergence of topological as well 
as dynamical nontrivialities is offered by Fig. [7] where the same separation 
of a strong-coupling piece of the domain is shown to survive up to the 
very extreme of c « 1. 

4 The domain of cryptohermiticity in the full- 
fledged three-parametric dynamical scenario 
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Figure 8: Boundaries Q(a, b, c) = and forbidden parts of the a — c plane as 
sampled at b 2 = 1. 



4.1 The auxiliary domains and their boundaries 

The domain T>^ of parameters for which our Hamiltonian H( 4 '(a, b, c) gen- 
erates the unitary evolution of quantum system is defined as an intersection 
of the triplet of domains T>^ p,( ^ ,w ^ in M 3 . After the preparatory considera- 
tions as presented in the preceding sections let us now leave all of the three 
parameters a, b and c freely variable and recall that 

• the domain T)^ is defined by the inequality P(a,b,c) = 10 — a 2 — 
2b 2 — c 2 > 0. It is compact so that we may restrict our attention just 
to the intervals of b 2 < 5, a 2 < 10 and c 2 < 10. At any fixed b 2 < 5 
the section of this first auxiliary domain coincides with the interior of 
a central circle in a — c plane with radius R = a/10 — 2b 2 ; 

• the allowed interior of the triply connected domain T>^> is defined 
by the inequality Q(a, b, c) = [(a + 3)(c - 1) - b 2 ] [(a - 3)(c + 1) - b 2 ]. 
In the a — c plane the boundaries of this domain are two hyperbolas 
sampled at b 2 = 1 in FigJHJ 

• the third auxiliary domain X)W is defined by the inequality W(a, b, c) = 
(8 + c 2 - a 2 ) 2 - 4 [16 - (a + c) 2 ] b 2 > 0. 
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The description of the latter domain D^' is slightly less trivial. It may 
be based on the observation that the interior of this domain covers all the 
exterior of the strip where \a + c\ > 4. Then the interior of this strip may be 
reparametrized, 

c — a = 2 r(c, a) G (—00, 00) , c+a = 4 sin <p(c, a) , ip(c, a) G (— ir/2, tt/2) 
making the rest of the domain T>( w > determined by the elementary inequality 

^ < \^ + r(c,a)Mv(c,a)]\ m (n) 
cos[(^(c, a)\ 

This means that within the restricted range of r(c, a) G (— vTo, vTo) the 
growth of |6| — > 00 must be compensated by the decrease of cos[<^(c, a)] — > 0, 
i.e., by the convergence c — > ±1 — a. This makes the strip- restricted part of 
the domain very small but increasing with the decrease of 1 6 1 from a 

sufficiently large initial value. 

4.2 The boundaries of the cryptohermiticity domain 

In the light of the latter comment it makes good sense to start the study 
of the "physical" overlaps of the three auxiliary domains J)( p >Q> w ) at 

the maximal admissible plane of b = = \/5 which touches the boundary 
at a = c = 0. This points still lies out of the domains P^' and T)^ 
since W(0, a/5, 0) = 64 — 320 < 0. In a search for the first touch between 
the b— plane and boundary we must diminish our b and move into the 

interior of T>^ p \ 









-3 
-2 






-0 c 

— 1 






—2 
-3 


-3-2-10 1 2 3 



Figure 9: b = VE - 1/100 = 2.226067977. 

In the first illustrative example at b = a/5 — 1/100 our Fig. [9] displays the 
motion of the triplet of boundaries dT>( p 'Q' w ^ projected into the a — c real 
plane. This picture shows that the corresponding section of the first domain 
becomes nonempty. Still, it just occupies the interior of a very small 
circle C(b) = dV^\ b= f ixed with the center at the origin. 

The interior of the second domain T>^> is perceivably bigger since, in 
the manner indicated by Fig. [8] above, it occupies the large domain between 
the two outermost, b — dependent hyperbolic curves Hi )2 (6) C dT>^ '\b= fixed- 
The overlap T)^ itself remains empty because the third domain V^ w ^ is 
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Figure 10: b = y/E - 1/2 = 1.736067977. 




Figure 11: b = V$ - 1 = 1.236067977. 

localized behind the two remaining and less trivially parametrized curves 
Gi, 2 (6) C dT>W\ b 

= fixed' 

During the subsequent decrease of b sampled by Fig. the two curves 
Gj(b) and H 3 -(6) (assigned the same subscript j — 1 or j — 2) get closer 
to each other while the internal circle C(b) gets larger. At each j and at 
the same value of b both the curves Gj, Hj touch the circle C(b). At a still 
smaller 6 = V5 - 1/2 = 1.736067977 they already move inside, sharing their 
two separate intersections with the circle. This situation is illustrated in 
Fig. EH 




Figure 12: b = 1.01. 
At any j and during the further decrease of the parameter b the two 
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Figure 13: b = 1. 



triple- intersection points between C(b), Gj and H., move apart. Between 
them one discovers the formation of the first two non-empty components of 
the physical domain T>^ H \ These components are disconnected, extremely 
narrow and eye-shaped, with the "eyes almost closed" but "slowly opening" 
with the further decrease of b. Graphically, the generic situation of this type 
in illustrated by Figs. [LTJ and \¥2[ 




Figure 14: b = 0.999. 

The next qualitative change of the pattern occurs at the above-mentioned 
special value of b = 1 as which we touch the saddle of the surface &D^ W >. 
Slightly before this happens we encounter the situation depicted in Fig. [T2l 
where the two separate subdo mains of the physical domain T>^\ b= f ixed al- 
ready almost touch. Next they do touch (cf. Fig. [T3]) and, subsequently, get 
connected (cf. the next Fig. [T41) . 

Surprisingly enough, below the saddle point 6=1 the topological sur- 
prises are still not at the end. There are no real news even at b = 0.6 (cf. 
Fig. [15]). Nevertheless, in the latter picture we already must pay attention 
to the two subdomains with the maximal a 2 s. 

Having selected just the right end of the (symmetric) picture at the pos- 
itive a ~ 3 we reveal the emergence of a tendency towards a new intersec- 
tion between the (up to now, safely external and non-interfering) second 
branches of the Q(a, b, c, )— related hyperbolas jj^ econd ' ) and of the back- 
bending boundaries dT>( w ^\ b= f ixed . For example, these curves get very close 
to each other but still do not intersect yet at b = 0.4, staying also outside of 
the central circular domain (cf. Fig. [T6|) . 
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Figure 15: b = 0.6. 




Figure 16: b = 0.4. 



The change of the pattern is finally achieved slightly below b = 0.4, at 
the moment when both of the new intersection candidates touch the circle 
C(6) = dV^\ b= f ixed in a single point. Subsequently, this point splits into 
the pair of the triple intersections and the further decrease of |6| forms the 
pattern which is sampled in Fig. ITTlat b = 0.2 and in Fig. HBlat b = 0.1. 

Certainly, the key features of the new situation are better visible, at the 
illustrative b = 0.1, in its magnified presentation as mediated by Figs. [[9] 
and l20l) . Now we may return to the limiting pattern of Fig. H] where we 
witness the abrupt change of the topology caused by the final confluence of 
the straight and backbanding branch of the boundary &D^ w '\f,= fixed in the 
limit 6^0. 




Figure 17: b = 0.2. 



12 



Figure 18: A return to the triply connected J)^ H ' at b = 0.1. 




2.6 2.7 2. a 2.9 3 3.1 3.2 



Figure 19: A magnified detail of Fig. [T8j 

Obviously, in a way confirmed by the complementary results of section 
13.31 such a confluence of the boundaries only occurs in the limit so that the 
thee-dimensional version of the open set T)( H > is ultimately confirmed to be 
triply connected. 

5 Summary 

In the history of pure mathematics the specification of the horizons dT>^ H > 
(called, often, "discriminant surfaces" in this context) has been perceived 
as a challenging and rather difficult problem even in its first "unsolvable" 
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Figure 20: A magnified detail of Fig. [19j 
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case characterized, in our present notation, by the Hilbert-space dimensions 
N = 5 It is also certainly worth noticing that in certain mathematically 
natural directions the real progress is of an amazingly recent date [12]. In 
this connection it seems rather remarkable that certain parallel developments 
occured also in several applied-mathematics oriented studies paying attention 
to the natural presence of more symmetries in the Hamiltonian \7\ and/or 
to the introduction of more observable quantities within a phenomenological 
quantum model in question [13]. 

In a constructive, more pragmatic setting as sampled by our recent paper 
[H] we restricted our attention to the topological aspects of the problem of 
horizons. The most obvious motivation of such an effort has been given by 
the fact that a disconnectedness of the domain T>( H ' immediately requires 
the transition from its traditional perturbation-theory descriptions (with a 
recommended recent compact sample given in [14J) to non-perturbative meth- 
ods, or to certain perturbative strong-coupling techniques at least [T5] . In 
Ref. [6] the parallel and less formal motivation has been emphasized to lie in 
the systematic search for possible physical origin of the dynamical anomalies 
in a kinematical nontriviality of the topology of phase space. 

The conclusions of our present paper are encouraging. Firstly we demon- 
strated that for many purposes it may be sufficient to use the matrices with 
a not too large N. Secondly, we have shown an increase of the feasibility 
of the description of models H with certain additional symmetries. At the 
first nontrivial Hilbert-space dimension N = 4 they reduced, for example, 
the minimal necessary number of parameters to D = 3 or even to D = 2 . 

Thirdly we clarified that once one works with a tridiagonal-matrix form 
of the "unperturbed" N hy N Hamiltonian Hq N ^ complemented by certain 
computationally suitable specific perturbations, the existence of the discon- 
nected subdomains in T>( H > opens the direct access to the strong-coupling 
dynamical regime. 

Fourthly, on mathematical side, we are now able to recommend, strongly, 
the use of certain auxiliary symmetries in the Hamiltonians (e.g., the ones 
recommended in Ref. [16J). In such cases, the algebraic secular equations 
pertaining to the model often happen to factorize, leading to the much better 
tractable polynomial equations of perceivably lower orders. 

Last but not least it seems worth emphasizing that even at our present 
"solvable" choice of the dimension N = 4 the latter fact rendered our toy 
model more easily solvable. In fact, on the background given by Refs. [SI ELT] 
it took some time for us to imagine that the anomalies of spectra could cer- 
tainly start occurring at the very small dimensions N < 8. In this context 
the message of our present N = 4 study is encouraging. Partly because sev- 
eral specific spectral irregularities as observed at N = 8 in Ref. [Bj do also 
exist for certain very similar Hamiltonian matrices with the dimension as 
low as iV = 4, and partly because our model reconfirmed the hypothesis of 
a "hidden" topology-related connection between the loop-shaping of the lat- 
tices (i.e., presumably, Betti numbers in continuous limit) and the existence 
of definite strong-coupling dynamical anomalies in the spectra of the energy 
levels. 
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Appendix A. The three-Hilbert-space formu- 
lation of quantum mechanics 

In section 2 of Ref. [5] one finds one of the most compact introductions into 
the abstract formalism of VT— symmetric quantum mechanics (PTSQM). 
Thus, we may shorten here the introductory discussion and refrain ourselves 
to a few key comments on the general theoretical framework. 

In such a compression the PTSQM formalism may be characterized as 
such a version of entirely standard quantum mechanics in which, in princi- 
ple, the system in question is defined in a certain prohibitively complicated 
physical Hilbert space of states where the superscript may be read as 
abbreviating "prohibited" as well as "physical" [3]. 

Typical illustrative realistic examples may be sought in the physics of 
heavy nuclei where the corresponding fermionic states are truly extremely 
complicated. In the latter exemplification the first half of the PTSQM recipe 
lies in the transition to a suitable, unitary equivalent Hilbert space, — > 
1-L^ S \ where the superscript may stand for "suitable" or "simpler" [3]. 

In the above-mentioned realistic-system illustration, for example, the new 
space coincided with a suitable "interacting boson model" (IBM). In the 
warmly recommended review paper of this field [T5] it has been emphasized 
that the requirement of the unitary equivalence between the two Hilbert 
spaces •H( p ) and may only be achieved in two ways. Either the cor- 
responding boson-fermion-like mapping Q between these two Hilbert spaces 
(known, in this context, as the Dyson's mapping) remains unitary (and the 
mathematical simplification of the problem remains inessential) or is admit- 
ted to be non-unitary (which is a less restrictive option which may enable 
us to achieve a really significant simplification, say, of the computational 
determination of the spectra). 

What remains for us to perform and explain now is the second half of 
the general PTSQM recipe. Its essence lies in the weakening of the most 
common unitarity requirement imposed upon the Dyson's mapping, 

fi 1 " = Or 1 

to the mere quasi-unitarity requirement 

fit = en -1 . 

The symbol 8^ / represents here the so called metric operator which defines 
the inner product in Hilbert space 

More details using the present notation may be found in [3]. Just a few of 
them have to be recalled here. Firstly, the main source of the purely technical 
simplifications of the efficient numerical calculations (say, of the spectra of 
energies) is to be seen in the introduction of the third, purely auxiliary Hilbert 
space where the superscript combines the meaning of "friendlier" 
with "falsified" [3]. 

By definition, the two Hilbert spaces and 'HS F > coincide as the mathe- 
matical vector spaces ("of ket vectors" in the Dirac's terminology). One only 
replaces the nontrivial metric 0^ = VftVL of the former space by its triv- 
ial simplification = I in the latter Hilbert space. As an immediate 
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consequence, the latter space acquires the status of an auxiliary, manifestly 
unphysical space which does not carry any immediate physical information or 
probabilistic interpretation of its trivial though, at the same time, maximally 
mathematically friendly inner products. 



Appendix B. The role of VT— symmetry 

In its most widely accepted final form described in Ref. [1] the PTSQM recipe 
complements the latter general scheme by another assumption. It may be 
given the mathematical form of the introduction of the second auxiliary, 
manifestly unphysical vector space KP 3 ^ which is, by definition, not even the 
Hilbert space. In fact, this fourth vectors space is assumed endowed with the 
formal structure of the so called Krein space [19] . 

Ref. [2] may be consulted for more details. Here, let us only remind 
the readers that the symbol V in the superscript carries the double meaning 
and combines the mathematical role of the indefinite metric V (defining in 
fact the Krein space) with an input physical interpretation (usually, of the 
operator of parity). In addition, the theoretical pattern 



(12) 



gets complemented by the requirement that there exists a "charge" operator 
C such that the (by assumption, non-trivial, sophisticated) metric 0^ 5 ^ ^ I 
which defines the inner product in the second Hilbert space H^' coincides 
with the product of the two above-mentioned operators, 



6 {s) = VC . 



(13) 



The contrast between the feasibility of the N ^ 3 constructions presented in 
Ref. [5] and the discouraging complexity and incompleteness of the next-step 
N = 4 constructions as performed in our older paper [17] and in its sequels 
[7J was also thoroughly discussed in our review [8j . In our present text we do 
not deviate from the notation and conventions accepted of this review. We 
solely pay attention to the class of models where N— dimensional matrix of 
parity is unique and given, in advance, in the following form, 



V 



' 1 









... 





-1 







... 








1 





... 











-1 












'•• 


. 










Tl 



(14) 



In parallel, we did make use of the time-reversal operator T of the form 
presented, e.g., in Ref. [5] as mediating just the transposition plus complex 
conjugation of vectors and/or matrices. One should add that once we solely 
work with real vectors and matrices, we are even allowed to perceive T as a 
mere transposition. 
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